February 2, 2008 
HU-EP-05-66 
DESY-05-211 

SFB/CPP-05-69 



Explicit results for the anomalous three point function 
and non-renormalization theorems 



F. Jegerlehner" and O.V. Tarasov'' 



" Humboldt Universitdt zu Berlin, Institut fiir Physik 
Newtonstrasse 15, D- 124-89 Berlin, Germany 
E-mail: f j egerSphysik . hu-berlin . de 
Deutsches Elektronen Synchrotron DESY 
Platanenallee 6, D- 15738 Zeuthen, Germany 
E-mail: 01eg.Tarasov@desy.de 



Abstract 

Two-loop corrections for the (VVA) correlator of the singlet axial and vector currents in QCD 
are calculated in the chiral limit for arbitrary momenta. Explicit calculations confirm the non- 
renormalization theorems derived recently by Vainshtein [Phys. Lett. B 569 (2003) 187] and Knecht 
et al. [JHEP 0403 (2004) 035]. We find that as in the one- loop case also at two loops the (VVA) 
correlator has only three independent form-factors instead of four. From the explicit results we 
observe that the two-loop correction to the correlator is equal to the one-loop result times the 
constant factor C2{R)as/Tr in the MS scheme. This holds for the full correlator, for the anomalous 
longitudinal as well as for the non-anomalous transversal amplitudes. The finite overall Us depen- 
dent constant has to be normalized away by renormalizing the axial current according to Witten's 
algebraic/geometrical constraint on the anomalous Ward identity [{VVdA) correlator]. Our obser- 
vations, together with known facts, suggest that in perturbation theory the (VVA) correlator is 
proportional to the one-loop term to all orders and that the non-renormalization theorem of the 
Adlcr-Bcll-Jackiw anomaly carries over to the full correlator. 



1 Introduction 



The Adler-Bell-Jackiw ^ |2] triangle anomaly in the divergence of the axial vector current is well 
known to play a crucial role at several places in elementary particle physics. Its nature is controlled by 
the Adler-Bardeen non-renormalization theorem 0, by the Wess-Zumino integrability condition and 
the Wess-Zumino effective action j4j, by Witten's algebraic/geometrical interpretation which requires 
the axial current to be normalized to an integer [S] and by the t' Hooft quark-hadron duality matching 
conditions [Sj. Phenomenologically, it plays a key role in the prediction of vr" — > 77, and in the solution 
of the U{1) problem. Last but not least, renormalizability of the electroweak Standard Model requires 
the anomaly cancellation which dictates the lepton-quark family structure. 

More recently Vainshtein ^7^ found an important new relation wxiq^) = ^Wbio^) (see below) matching 
to all orders in perturbation theory, in some kinematical limit, the transversal part to the anomalous 
longitudinal amplitude which is subject to the Adler-Bardeen non-renormalization theorem. Later 
Knecht et al. |F were able to generalize this kind of non-renormalization theorems. These recent 
investigations came up in connection with problems in calculating the leading hadronic effects in the 
electroweak two-loop contributions to the muon anomalous magnetic moment 1101 1111 112j . 
The first electroweak two-loop calculations 13 for revealed that triangle fermion-loops may give 
rise to unexpectedly large radiative corrections. The diagrams which yield the leading corrections are 
those including a VVA triangular fermion-loop (y^^ 7^ while VVV = ) associated with a Z 
boson exchange 




and a fermion of flavor / gives a contribution, up to UV singular terms which will cancel. 
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where a is the fine structure constant, the Fermi constant, T^f the 3rd component of the weak 
isospin, Qf the charge and Ncf the color factor, 1 for leptons, 3 for quarks. The mass nif is if 
ruf < rUfj, and ruf if mj > m^. Cf denotes constant terms. Since, as granted in the Standard Model 
of elementary particles, anomaly cancellation by lepton-quark duality J2f ^cfQ'jT^f = is at work, 
only the sums over complete lepton-quark families yield meaningful results relevant to physics. In 
any case the quark contributions have to be taken into account. In fact treating the quarks like free 
fermions (quark parton model QPM) the first family yields 
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which demonstrates that the leading large logs ~ In Mz have dropped! However, the quark masses 
which appear here are illdefined constituent quark masses, which can hardly account reliably for the 
strong interaction effects. 

Since we are interested in a static low energy quantity = ^ (5 — 2)^ = Fm(0), given by the Pauli 
form factor at zero momentum transfer, the perturbative QCD (pQCD) calculation of the light quark 
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contributions seems more than questionable. However, there is a large scale in the game, namely the Z 
boson mass Mz ~ 91.19 GeV, which makes an estimate by the quark parton model as a first step not 
completely unreasonable. Indeed, in the relevant kinematical region, the leading strong interaction 
effects may be parametrized by two VVA amplitudes, a longitudinal and a transversal wr(Q^) 

one, which contribute as 1111 1121 



VVA 
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where A is a cutoff to be taken to oo at the end. For a perturbative fermion loop to leading order [14j 

dx x(\ — x) 
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Vainshtein [7j has shown that in the chiral limit the relation 
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is valid actually to all orders of perturbative QCD. Thus the non-renormalization theorem valid beyond 
pQCD for the anomalous amplitude wl (considering the quarks q = u,d, s, c, b, t only): 
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(1.4) 



carries over to the perturbative part of the transversal amplitude. Thus in the chiral limit the per- 
turbative QPM result for wt is exact. This may be somewhat puzzling, since in low energy effective 
QCD, which encodes the non-perturbative strong interaction effects, this kind of term seems to be 
absent. 

We observe that the contributions from wl for individual fermions is logarithmically divergent, but 
it completely drops for a complete family due to the vanishing anomaly cancellation coefficient. The 
contribution from wt is convergent for individual fermions due to the damping by the Z propagator. 
In fact it is the leading l/Q^ term of the wt amplitude which produces the In^^ terms. However, 
the coefficient is the same as for the anomalous term and thus for each complete family also the In Mz 
terms must drop out. 

Low energy QCD is characterized in the chiral limit of massless light quarks u,d,s, by spontaneous 
chiral symmetry breaking (S^SB) of the chiral group SU{3)v SU{3)a, which in particular implies 
the existence of the pseudoscalar octet of pions and kaons as Goldstone bosons. The light quark con- 
densates are essential features in this situation and lead to non-perturbative effects completely absent 
in a perturbative approach. Thus such low energy QCD effects are intrinsically non-perturbative and 
controlled by chiral perturbation theory (^PT), the systematic QCD low energy expansion, which 
accounts for the SxSB and the chiral symmetry breaking by quark masses in a systematic manner. 
The low energy effective theory evaluation of the hadronic contributions related to the light quarks 
u, d, s was worked out in [21^3 and later in |J^. It was shown that in the operator product expansion 
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(OPE) the leading non-perturbative (NP) term in the chiral Hmit is due to the u, d, s quark condensate 



^.r(g2)^P^-vr2_-il^ at Q2 large 
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which breaks the degeneracy WTiQ"^) = \ wiiQ"^) found in perturbation theory^. Mp is the p mass. 
Relevant effective couplings are the neutral current part 
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2 sin Qw cos Qw 
and the Wess-Zumino Lagrangian 



A9^(7r0 + i=r?8-i=%) Z^ 
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which reproduces the ABJ anomaly via the PCAC relation, e is the positron charge, sin^ is the 
weak mixing parameter, the pion decay constant and vrO is the neural pion field. The pseudoscalars 
r/8,r/o are mixing to ?y, The [n, d, s] contribution evaluated this way is given by the diagram 




which together with its crossed version in the unitary gauge and in the chiral limit and completed for 
the first two fermion families yields jHl 
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The residual M^-dependence in effective theory was controversial^, and is contradicting Vainshtein's 
non-renormalization theorem [71 In fact, there is not really a contradiction. The point is that 
the low energy effective theory does not apply up to the Mz mass scale as assumed in obtaining the 
above result. One rather has to apply a cut-off Amatching of the order of the proton mass m^, say, and 
above that scale of course the QPM gives the correct answer. The leading large log terms proportional 
to In Mz then again cancel and the M^-dependence gets replaced by a Ainatching~dependence (with 
coefficient -14/3 changed to -2/3 and -35/3 to -47/3) in the above result. 

As the extensions of the Adler-Bardeen non-renormalization theorem for the anomalous Ward identity 
(VVdA) turn out to play an important role in new phenomenological applications, we will study in 
the following such possible generalizations by an explicit calculation of the leading QCD corrections 
to the j'yZ triangle. 

^The OPE only provides information on wt for large. At low we only know that 10^(0) = 1287r^ where 
C22 is one of the unknown xPT constants in the 0{p^) parity odd part of the chiral Lagrangian 1151 . The low energy 
effective theory does not yield a 1/Q^-pole as required by the non-renormalization theorem. In the OPE terms 
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are due to explicite chiral symmetry breaking which yields AwriQ )np 

■^Numerical estimates presented in and |12) agree well within errors and the discrepancy in question is of conceptual 
nature and not relevant for the interpretation of the present g — 2 experiment. 
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2 Definitions 



Let us consider the VVA three point function 

W^,p(gi,(?2) = i / d*xid^X2e'^'i^-^'+i^-^^^ X {0\T{V^{x^)V,{x2)ApiO)} \0) (2.6) 



of the flavor and color diagonal fermion currents 
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where tp is a quark field. The vector currents are strictly conserved d^jV^ = 0, while the axial vector 
current satisfies a PCAC relation plus the anomaly df^A^^ = 2imoV'75V' + ^^fj.upaF^'^ {x)FP" {x). We 
will be mainly interested in the properties of strongly interacting quark flavor currents in perturbative 
QCD. To leading order the correlator of interest is associated with the one-loop triangle diagram 
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plus its crossed ((?i,yu Q2,t^) partner. In the following we will closely follow the notation of "HI. 
The Ward identities restrict the general covariant decomposition of yVf_iupiQi,Q2) into invariant func- 
tions to four terms 

+ ^5^^ (g?,gi,gi) t^ptliqi,Q2) + w'f'' (ql,ql,qt) tl'Jpiqi,q2) + {ql,qi,ql) iti(9i>92)| ,(2.8) 

with the transverse tensors given by 
*K(gi,92) 
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Bose symmetry ((71,/x <-> q2-,v) entails 
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The longitudinal part is entirely fixed by the anomaly, 

WL (qhqhqt) = 



2N, 
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which is exact to all orders of perturbation theory, the famous Adler-Bardeen non-renormalization 
theorem. 

In (5j the following three chiral symmetry relations between amplitudes were derived in pQCD: 
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involving the transverse part of the {VV A) correlator W^i/p(gi, ^2)) and which hold for all values of 
the momentum transfers q\, g| and g|. In the kinematical configuration relevant for g-2 calculations, 
gi = /c ± g, q2 = —k, expanding to linear order in k and noting that i^ti/p((?i, 52) ~ ^Attp(9i) 92) = 
q^eiiupak'' - qf_ieypapq°'k^ - qpS^ua/sq'^k'^ + 0(A;^) and t'i^l{qi,q2) = 0{k'^), these relations imply the 
non-renormalization theorem Hl.cil) obtained in Refs. upon identifying 

wl{Q^)=wl{-Q^O,-Q^) , WTiQ^)=w^Tt\-Q^,0,-Q^) + w^-\-Q\0,-Q^) , (2.13) 
with Q'^ = —q^. 

3 Calculations 

We perform the calculation with conventional dimensional regularization 'IT and use a linear covariant 
gauge with arbitrary gauge parameter ^ throughout the calculation. 



A 




V„ 









Two-loop QCD diagrams contributing to (yVA) correlator 
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Our procedure of treating 75 is similar to the one used in jl6j . We write down all fermion loops 
starting with the axial-vector vertex, and then perform Feynman integrals and Dirac algebra without 
assuming any property of 75 at all. In this way all diagrams will be expressed in terms of traces of 10 
combinations of 7 matrices: 



4iAi = 7p757m7i^^2, 4i^2 = IplbQiliilu, "iiAs = ^2 7^75^17/^92, 

4iA4 = Qi7p75gi7/^92, 4i^5 = -q^ilpimHi^h, ^^i^e = -Qilplbqiltih, .3 

HAj = 'jpj^qi, 4z^8=7p75Q2, iiAg = 'jp^^'j^, 

AiAio = jp-f5ju 

with q = qpj^. The prescription is sufficient to enable us to arrive at expressions in front of ^di, . . .Aiq 
which have finite limits as d ^ 4. After this the usual formulae 

Tr [75 7„ 7/37^7,,] = 'iieafipu , Tvlj^-falp] = (3.15) 

valid in d = 4 dimensions were used. In our convention £0123 = +1 and (1 — 75)/2 projects to 
left-handed fermion fields. 

Tensor integrals were expressed in terms of integrals with different shifts of the space-time dimen- 
sion IJH] . All scalar integrals were reduced to 6 master integrals by using the Grobner basis technique 
proposed in ^Hl- The expressions for the individual diagrams are sums over 15 terms which are 
combinations of the 6 basis integrals 
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kl{ki-qir{ki-q2r 
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A'Hq]) = 

R2{q!,qlqh = 

P.{<lUl4) = y y P;2p\2fc2(^^_..^2^: ..2n.^ , ._^2 ■ (3-16) 
and multiplied by ratios of polynomials in momenta and d: 

^1 Qk[qi,a) 

The momentum dependence of the denominators turns out to be rather simple: 

Qkiqld) = Q{d){qlr{qlf^{qlr/S.''^ (3.18) 
where ak,bk,Ck,ek, are some numbers, Q{d) is a polynomial in d and 

A = qf + ql + ql- Iqjqj - 2qlql - 2qjql . (3.19) 
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The integrals H3.16|) form a complete set of master integrals needed for the calculation of massless 
vertex diagrams with planar topology. 

The sum of all diagrams turns out to be gauge parameter independent. 

In the Feynman gauge at (73 = and for arbitrary d, the results of our calculation are in agreement 
with the ones presented in [l^ diagram by diagram. 

By applying the prescription outlined above for the evaluation of individual diagrams, the tensor 
structures 

Ni = q^EafS^iu , N2 = q2£a^3^lu , (3.20) 

which appear in the Ansatz (|2.8j) are actually absent. In place of Ni,N2, two other tensor structures 
are present which we eliminate by means of the Schouten identities: 

^ll2^1^a(3yp = —QlQ2Qip^al3pu + ^2 9l + (ll(ll(l2^apvp — (l2(ll^ap.up ■ (3.21) 

Reshuffling terms in this way allows us to express each diagram in terms of the tensor structures 
introduced in (|2.8)) exhibiting manifestly the vector current conservation. 



4 Results and Discussion 

Including one- and two-loop contributions, we may represent the form-factors in the form 

w'r^iql^ = rif NcW^^^{ql,q2,ql) + a nj NcC2{R) w^^r^{ql,q2,ql) 

WLiqlqlqh = N^wi^L{ql,ql,ql) + a Uf N^C2{R) W2,L{ql,ql,ql) (4-22) 



where 



includes the QCD coupling a^, g as usual is the gauge coupling, n/ is the number of flavors and Nc the 
number of colors. The quarks are in the fundamental representation R and the corresponding group 
theory factor is given by 

C2{R)I = R^R" , C2(i?)=4/3 for QCD . (4.24) 

We have been working in the MS renormahzation scheme. The singlet axial current = Ap is non- 
trivially renormalized because of the axial anomaly. It is known fl^ that in addition to the standard 
ultraviolet renormahzation constant which reads = 1 in our case (as Zj^ — 1 = 0(a^)), one 
has to apply a finite renormahzation constant Z^ such that renormalized and bare currents are related 
as: 

{j'p)r = Z,Zj^{j'p)o. (4.25) 

The counterterms coming from the wave function renormahzation of quarks and ultraviolet renormal- 
ization of the axial and vector currents cancel. The finite renormahzation constant is known at the 
three- loop level -^o^ calculations we take 

Z5 = l-4C2(i?)a . (4.26) 
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The result of the two-loop calculation after adding all diagrams is surprisingly simple and, normalized 
according to H4.22() . is given by 



qiw2,L{<lhQ2,li) = -8 (4.27) 

4~Ti<ll^<llQl) = -w^^kllllll), (4.28) 

ql/:^^w':^kqlqlql) = 8(x - 2/)A + 8(x - y)(6xy + A)$« (x, y) 

- 4[18x?/ + Qx^ -Qx + {l + x + y)A)]La; 

+ 4[18xy + - 6y + (1 + X + y)l^)]Ly (4.29) 

ql^^w^^^qlqlql) = 8[6xy + (x + y)A]$« (x, y) + 8A 

- 4[6x + A](x-y-l)L^ 

+ 4[6y + A](x-y + l)Ly (4.30) 



with 



Lx = \nx, Ly=lny, x = V = ■ (4-31) 

^3 93 

The explicit expression for may be found in j25j : 

cl>(i)(x, y) = ^ I 2 (Li2 (-px) + Li2 i-py)) + In I In + In(px) In(py) + ^ | , (4.32) 

where 

A(x,?/) = VA , p(x,y) = 2 (l-x-y + A)-\ A = (1 - x - y)^ - 4xy . (4.33) 
The comparison with the results of the one-loop calculation reveals that 

y^lipp{(ll,(l2) \two-loop = 4.C2{R)aW^up{qi,q2) \one-loop (4.34) 

Multiplying the sum of one- and two-loop terms by the finite factor Z5 we arrive at 

) 92) I one— loop 

(4.35) 

This is the non-renormalization theorem for the full off shell correlator at two-loops. While our 
calculation confirms the relations (|2.12|) derived in 8_ and the non-renormalization theorem H1.3|) 
found in j3 El , these findings are not sufficient to explain our result valid for generic momenta. 
Taking into account the rather non-trivial momentum dependence of the form-factors it is very tempt- 
ing to suggest that it could hold to all orders of perturbation theory because of the topological nature 
of the anomaly, for example. 

We would like to stress that the surprising relation could be discovered only by keeping the gen- 
eral non-trivial momentum dependence. The anomalous three point correlator exhibits an unusually 
simple structure, while contributions from individual diagrams are very unwieldy. One can expect 
similar effects for other anomalous correlators. Since at the order considered the QCD calculation is 
essentially a QED calculation, it is highly non-trivial whether this carries over to higher orders. For 
the (VVdA) anomalous correlator a large number of two-loop calculations have been performed ((26j- 
j32j). mainly in QED and we refer to the comprehensive review by Adler |33j and the references therein. 
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In electroweak SM calculations one would a priori expect that remormalizing parameters and fields 
would be sufficient for renormalizing the SM. Our calculation shows that on top of the standard renor- 
malization, it is mandatory to renormalize the anomalous currents by the finite renormalization 
factor Z5 because the lepton currents and the quark currents pick different Z-factors and if they 
are not renormalized away the anomaly cancellation and hence renormalizability obviously would get 
spoiled. As pointed out by Adler and many others j33j the point is there exists a renormalization 
scheme for which the one-loop anomaly is exact. Only in this scheme anomaly cancellation and thus 
renormalizability will carry over to higher orders in the SM. Our result shows that due to the necessity 
of renormalizing away possible higher order contributions from the anomaly also the non-anomalous 
transversal contributions are affected. We have shown that at least at two-loops the entire contribu- 
tion gets renormalized away in the zero mass limit. 
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